BE IT KNOWN We. Tobias LANG and Uwe KONZELMANN, 

have invented certain new and useful improvements in 

METHOD FOR PRODUCING AT LEAST ONE CHARACTERISTIC LINE 
OF AN AIR MASS DETECTING DEVICE FOR AN INTERNAL 
COMBUSTION ENGINE 

of which the following is a complete specification: 
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BACKGROUND OF THE INVENTION 



The present invention relates to a method for producing at least 
one characteristic line of an air mass detecting device for an internal 
combustion engine. 

Such methods are known in the art. In this method an air mass 
detecting device on the one hand and an accurate comparison probe on the 
other hand are installed in an impeller testing stand which must simulate the 
flow situation in an aspiration region of an internal combustion engine. The 
signals of both sensors are extracted. From the output signals of the air 
mass detecting device and from the mass stream detected by the 
comparison probe, a characteristic line is formed which takes into 
consideration the influence of the geometry of the aspiration region of the 
intemal combustion engine and its action on the signal. 

In normal operation of internal combustion engine the signal of 
the air mass detecting device serves, in addition to other criteria, for 
determination of the loading condition of the internal combustion engine. 
Conventionally, a hot film air quantity sensor is used as the air mass 
detecting device, which is Identified as "HFM-sensor". 



The characteristic line determined by the testing stand research 
is identified as a "static" characteristic line, since it is produced at static or 
stationary flow conditions. It is provided in a control device of the internal 
combustion engine. The problem is however that in many internal 
combustion engines the airflow in the aspiration region is not stationary, but 
instead has a pulsed pattern. The correct detection of such a pulsing airflow 
in real use is difficult for conventional air mass detecting devices in principle, 
so that a faulty indication is produced which is a function of the frequency 
and the amplitude of the flow pulsations. 

When however the air mass screen determined by the air mass 
detecting device does not correspond to the actual air mass stream in the 
combustion chamber of an internal combustion engine, deviations of the 
emission behavior of the internal combustion engine from an optimal 
emission behavior can take place. 



SUMMARY OF THE INVENTION 



Accordingly, it is an object of the present invention to provide 
a metliod of the above mentioned type, which provides for the air mass 
detecting device in condition of pulsating air flow such results which are as 
precise as possible. 

In keeping with these objects and with others which will 
become apparent hereinafter, one feature of the present invention resides, 
briefly stated, in a method of the above mentioned type, in which a dynamic 
characteristic line is produced in a method which includes the following 
steps: 

b) receiving signals of an air mass detecting device, by placing the air 
mass detecting device on an Internal combustion engine testing stand 
in different operational points of an air mass flow and detecting 
signals produced by the air mass detecting device, 

c) converting the signals of the air mass detecting device Into air mass 
flow values by interpolation on an output characteristic line; 

d) forming average values of the air mass flow values via integral 
plurality of a pulsation period for the corresponding operational point; 



e) calculating a deviation, which corresponds to the deviation of the 
average value of the air mass flow from an accurate comparison air 
mass flow, for the corresponding operational points; 

f) calculating the square norm through the matrix of the deviation; 

g) producing an adjusted characteristic line in the sense of an 
optimization with respect to the condition that the square norm is 
minimal. 

h) recalculating the signals of the air mass detecting device into air mass 
flow values by interpolation of the adjusted characteristic line; and 

i) iteration by repeating the steps h), d), e), f), g). 

The above mentioned objective is also achieved in a method 
of the above mentioned general type, in which a dynamic characteristic line 
is produced by a method which includes the following steps: 

a) obtaining the signals of the air mass detecting device, by placing the 
air mass detecting device in different operational points on an internal 
combustion engine testing stand of an air mass flow, and detecting 
signals produced by the air mass detecting device; 
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b) producing a histogram from the signals through at least one complete 
pulsation period for the corresponding operational points of the 
internal combustion engine; 

c) recalculating equidistant signal values into air mass flow values by 
interpolation on an output characteristic line; 

d) forming weighted average values of the air mass flow by using the 
histograms, correspondingly for the operational points; 

e) calculating a deviation which corresponds to the deviation of the 
average air mass flow from a comparison air mass flow 
correspondingly for the operational points; 

f) calculating the square norm through the matrix of the deviation; 

g) producing an adjusted characteristic line in the sense of optimization 
with respect to the condition that the square norm is minimal; 

h) recalculating the signals of the air mass detecting device into air mass 
flow values by interpolation of the adjusted characteristic line; and 

i: iterating by repeating the steps h), c), d). e), f), and g). 

The adjusted dynamic characteristic line formed in accordance 
with the inventive method leads, in particular in condition of significantly 
pulsating air flow, to higher accuracy in the determination of the actual air 
mass flow reaching a combustion chamber of an internal combustion engine 



from the output signal of the air mass detecting device. Finally, with this 
method the consumption and emission behavior of the internal combustion 
engine are significantly improved, since the mixture control is possible with 
a higher precision. 

With the adjusted dynamic characteristic line which deviates 
from conventional static characteristic lines, the dynamic flow behavior in an 
aspiration region of an internal combustion engine can be considered very 
well. In addition, corresponding data in a real internal combustion engine 
testing stand can be obtained in different operational points and stored. For 
example rotary speed and load, wherein the load can be indicated for 
example by the torque or the combustion chamber average pressure. With 
obtaining the data by a suitable selection of the operational points, 
advantageously the whole operation region of the internal combustion engine 
can be covered. 

The modified characteristic line in addition is independent from 
the presence of correction characteristic fields in a control device, which 
processes the output signals of the air mass detecting device, so that the 
inventive method can be used also with such control devices which do not 
have such a correction characteristic field. 



For forming an average value through the standard deviations 
in the testing stand, the signal of the air mass detecting device is indicated 
in the same time window as the signal of the comparison probe. This high 
time resolution of the measurement must not be worsened, since otherwise 
the important dynamic effects no longer can be correctly detected. For this 
reason when the inventive method a great grade quantity of data is 
developed. 

With the use of histograms proposed in the second-mentioned 
inventive method, a significant reduction of the data quantity with maintaining 
of the dynamically relevant informations is provided. For the corresponding 
operational point, instead of the full signal of the air mass detecting device, 
only a histogram of the signal is supplied. During the characteristic line 
optimization, it is not necessary to recalculate each individual signal value of 
the air mass detecting device by interpolation of the characteristic line in an 
air mass flow, but instead the limits of the equidistant histogram channels (in 
the practice its unit is volt) are replaced by a Binning- vector (unit for example 
kg/h) interpolated on the characteristic line. The decisive feature is that only 
a fixed number of interpolation is required, which is equal to the dimension 
of the Binning-vector. The number of the interpolations is also independent 
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from the scope of the measuring data, while the measuring accuracy 
improves with the scope of the measuring data. 

With the use of histograms it is possible to reduce the data 
quantity by several orders. This leads to the situation that the optimization 
algorithm converges significantly faster. Depending on the characteristic of 
the aspiration region of the internal combustion engine, with this method a 
convergence of an optimization process is possible. 

In accordance with the first embodiment of the invention it is 
proposed that the Levenberg-Marquardt method is used for the non-linear 
optimization. The non-linear optimization method converges relatively fast 
and is simple to program. Alternatively, genetic algorithms or evolution 
strategies for optimization can be utilized. 

The iteration can be intermpted after a predetermined number 
of iteration steps. Thereby the calculation expense remains within a 
predetermined range. 



Alternatively, it is also possible to interrupt the iteration at 
reaching a predetermined value for the square norm. In this case the 
accuracy of the optimization results is provided. 

It is further advantageous when different accidently generated 
static characteristic lines are used as output characteristic lines. This makes 
possible the recognition of sub-optimal extrema. The result of the 
optimization is thereby again improved. 

Further it is proposed to perform the optimization also with 
respect to a secondary condition, through which a desired course of the 
adjusted characteristic line is taken into consideration. Thereby for example 
a demand in accordance with a monotonous course of the adjusted 
characteristic line can be provided in calculations. 

The novel features which are considered as characteristic for 
the present invention are set forth in particular in the appended claims. The 
invention itself, however, both as to its construction and its method of 
operation, together with additional objects and advantages thereof, will be 
best understood from the following description of specific embodiments when 
read in connection with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a view schematically showing an internal 
combustion engine with an air mass detecting device; 

Figure 2 is a characteristic line of the air mass detecting device 

of Figure 1 ; 

Figure 3 is a view schematically showing a first embodiment of 
a method for producing a modified dynamic characteristic line for the air 
mass detecting device of Figure 1 ; 

Figure 4 Is a view showing two diagrams which provide data 
reduction by formation of histograms; 

Figure 5 Is a view showing a diagram, in which surfaces with 
the same relative deviation of the air mass flow determined from the statistic 
characteristic line from actual air mass flow are plotted at different 
operational points of the internal combustion engine of Figure 1 ; 



Figure 6 is a view showing a diagram which Is substantially 
similar to the diagram of Figure 5 on the basis of the adjusted dynamic 
characteristic line; and 

Figure 7 is a schematic similar to the view of Figure 3 for a 
second embodiment of a method for producing an adjusted dynamic 
characteristic line for the air mass detecting device of Figure 1 . 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 



A diesel internal combustion engine is identified in Figure 1 as 
a wliole witfi reference numeral 1 0. It includes several cylinders, from which 
only one is shown in Figure 1 for clarity. It has a combustion chamber 12, 
into which air is supplied through an aspiration pipe 14 and an inlet valve 16. 
The combustion gasses escape from the combustion chamber 12 through 
an outlet valve 18 and an exhaust gas pipe 20. The air mass flow which 
flows through the aspiration pipe 14 is detected by an air mass detecting 
device 24, which is formed in this case as an HFM sensor. 

The air mass stream is further detected by a high accuracy 
comparison probe 26, which is formed as a so-called "air clock" arranged in 
the exhaust gas pipe 20. Fuel is supplied directly to the Intemal combustion 
chamber 12 through an injector 28, which is supplied from a high pressure 
fuel system 30. An incandescent device 32 facilitates firing of the mixture in 
the combustion chamber 12 in condition of a cold start. 

The determination of the air mass which is supplied through the 
aspiration pipe 14 into the internal combustion chamber is very important for 
the correct mixture control in the combustion chamber 12. It is therefore 
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desired to detect the air mass flow by the HFM sensor 24 with highest 
possible precision. For this purpose a characteristic line is utilized which 
links the output signal Uhfm oiihe HFM sensor 24 with a corresponding air 
mass stream m. An example for such a characteristic line 38 is shown in 
Figure 2. It includes a plurality of support points 36. The characteristic line 
38 is produced by interpolation between the support points 36. Depending 
on the construction, more or less strong air pulsations can occur in the 
aspiration pipe 14. 

As a result of thermodynamic and aerodynamic effects on the 
HFM sensor 24, these air pulsations can lead to faulty measuring results, 
which can not be considered in the conventional use statistic characteristic 
lines. In order to minimize this error, a modified dynamic characteristic line 
is produced, which takes into consideration also dynamic flow effects in the 
aspiration pipe 14 and so exactly represent the air mass flow which flows 
through the aspiration pipe 14 to the combustion chamber 12. For this 
purpose a non-linear optimization process is performed, which is illustrated 
in Figure 3. 

First of all during a testing stand operation with the internal 
combustion engine 10, the coarse signals of the HFM sensor 24 are 
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indicated at different rotary speed/load points. In the shown example these 
signals are detected for 1 5 different rotary speeds and 1 5 different loads 60 
seconds long with a time resolution of 0.5 msec. This produces the output 
voltage Uhfm of the HFM sensor 24 as an array with the dimensions 15x15 
X 120000 (reference numeral 40 in Figure 3). 

This data quantity is reduced with maintaining the dynamical 
relevant information by the determination of histograms. For each rotary 
speed/load point of the Indicated time-dependent voltage signal U=f(t)(the 
uppermost diagram in Figure 4) for a complete pulsation period in a 
histogram n,erf(U)(central diagram in Figure4). A histogram corresponds on 
the one hand to an averaging over the whole measuring time t, or in other 
words over all standard deviations. On the other hand, it contains in it the 
total relevant dynamic information. 

The voltage Uhfm is represented in equidistant steps with a 
fixed step width (a region from zero to 5 volt is covered with a step width of 
0.005 volt).' The produced array n^, (reference numeral 42 in Figure 3) has 
in this example the dimension 15 x 15 x 1000. Thereby a reduction of the 
data quantity by two orders is achieved with approximately 180 MB to only 
approximately 1.8 MB. 
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In 54 the voltage values Uhfm of 0 to 5 volt (step width 0.005 
volt) is interpolated with square interpolation on a characteristic line provided 
in 52b. It has conventionally initially (Initial guess", reference numeral 51) 
distances A (reference numeral 52a) between the supporting points, 
which correspond to the distances between the supporting points (m, ) of 
a conventional static characteristic line. The above presented interpolation 
means nothing else but not each individual signal value of the HFM sensor 
24 is recalculated by interpolation on the characteristic line into an air mass 
flow, but instead however the limits of the equidistant histogram channels 
(unit'.volt) is replaced by a Binning-vector (unit:kg/h) interpolated on the 
characteristic line. The result is an air mass flow vector m hfm (reference 
numeral 56). 

In 44 for each operational point n, P^e a weighed average value 
for the air mass flow mHPM formed. It corresponds to the gravity point of 
the corresponding histogram. This leads to an air mass AowtThfm as a two- 
dimensional array, depending on the rotary speed n and the load P 
(reference numeral 46). In 48 from it the relative deviation dm/m is 
calculated which is the difference between the determined air mass flow 
(block 46) and the air mass flow mvs determined via the comparison probe 
26 referred to the air mass flow detected by the comparison probe (the air 
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mass flow is provided depending on the rotary speed n and the load P^e 
in 49 as a two-dimensional array). It is to be understood that in other 
embodiments, instead of the values determined by an accurate comparison 
probe, also the values determined in different ways can be utilized, which 
correspond the best to the actual air mass flow. 

With these matrix of the relative deviation dm/m, the square 
norm is now calculated in 50. It corresponds to the sum of the square of 
the deviations over all rotary speed-/load points. The calculation of the 
square norm is performed by a sums formation over the square matrix 
components. The optimization target is the minimization of this number, 
which leads to new distances A (reference numeral 52a) between the 
support points of the characteristic line. Therefore in 52b a modified 
characteristic line is obtained, which is characterized by the corresponding 
new supporting points m(air mass flow) and (voltage). 

Through the term X'^, secondary conditions can be taken into 
consideration during the optimization. For example a monotonous course 
can be extended from the characteristic line. Characteristic lines with non- 
monotonous course can be excluded as optimization result. This can be 
taken into consideration by a term X'^ which is great with a negative A U^. 
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By means of a square interpolation, this new characteristic line 52 is again 
recalculated to the equidistant voltages U^fm new supporting points which 
serve as Binning-vector and during the formation of the histograms. Thereby 
new air masses mHPM (reference numeral 56) are produced in accordance 
with the characteristic line in association with the voltages of the Binning 
vector. 

The steps 54, 56, 44, 46, 48, 50, 52a, and 52b are repeated in 
the sense of an iteration so often, until either a predetermined number of 
iteration steps is reached or the square norm X ^ reaches a predetermined 
value. Because of the data reduction by means of histogram, the time period 
required for the optimization amounts, in the presented example on a 
conventional calculation device, to approximately only 30 seconds. The time 
advantage produced by the data reduction is greater with the greater 
available data quantity. 

The accuracy of the modified dynamic characteristic line which 
is finally obtained in 52, when compared with a statistic characteristic line, is 
clear from a comparison of the diagrams of Figures 5 and 6. In them the 
surfaces are plotted with the same relative deviations dm/m, or in other 
words, the air mass flow determined by means of the corresponding 
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characteristic line from the actual air mass flow, depending on the rotary 
speed n and the average pressure P^g corresponding to the load. It can be 
seen that in the total operational region of the internal combustion engine 10, 
with the use of the modified dynamic characteristic line, the maximum 
relative deviation is between 6 and 10%, in wider regions is however only 
between -2 and +2% (Figure 6). With the use of a conventional static 
characteristic line, to the contrary at lower rotary speeds deviations up to 
18% are determined (Figure 5). 

Figure 7 shows a flow diagram of an alternative embodiment 
of the above described method. The same reference numerals as in Figure 
3 are utilized for functionally equivalent regions. The regions described in 
Figure 3 are not illustrated in connection with Figure 7 in detail. 

The method shown in Figure 7 differs from the method shown 
in Figure 3 in that a data reduction by formation of histograms is dispensed 
with. This leads to the situation that in 54 instead of 1 ,000 interpolations per 
rotary speed-/load point 27. millions interpolations (15 x 15 x 120000) are 
required. The corresponding calculation time for the optimization is therefore 
significantly longer than in the method shown in Figure 3 (approximately 6 
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weeks on the conventional computing device), and the optimization 
converges slower. 

It will be understood that each of the elements described 
above, or two or more together, may also find a useful application in other 
types of methods and constructions differing from the types described above. 

While the invention has been illustrated and described as 
embodied in method for producing at least one characteristic line of an air 
mass detecting device for an internal combustion engine, it is not intended 
to be limited to the details shown, since various modifications and structural 
changes may be made without departing in any way from the spirit of the 
present invention. 

Without further analysis, the foregoing will so fully reveal the 
gist of the present invention that others can, by applying current knowledge, 
readily adapt it for various applications without omitting features that, from 
the standpoint of prior art, fairly constitute essential characteristics of the 
generic or specific aspects of this invention. 
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What is claimed as new and desired to be protected by Letters 
Patent is set forth in the appended claims. 
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